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Summary 


An improved version of Woodward's chord plane aerodynamic panel 
method for subsonic and supersonic flow has been developed for cambered 
wings exhibiting edge-separated vortex flow, including those with leading- 
edge vortex flaps. The exact relation between leading-edge thrust and 
suction force in potential flow is derived. Instead of assuming the 
rotated suction force to be normal to wing surface at the, leading edge, 
new orientation for the rotated suction force is determined through con- 
sideration of the momentum principle. The supersonic suction analogy 
method is improved by using an effective angle of attack defined through 
a semi-empirical method. Comparisons of predicted results with available 


data in subsonic and supersonic flow are presented. 


1 . INTRODUCTION 


In references 1 and 2, an improved panel method was shown to be 
capable of predicting accurately the leading-edge and side-edge suction 
forces. In the method, a specific set of control point locations is 
obtained, based on a two-dimensional theory. All three-dimensional results 
presented for this method have been for non-carobered wings exhibiting 
edge-separated vortex flow. For highly swept cambered wings In subsonic 
compressible flow, a simplified method (as compared with that to be developed 
in this report) has been developed based on reference 3. That method 
uses the vortex- lattice method and suction analogy (VLH-SA) and is 
applicable only to subsonic flow. Its application to analysis and design 
of slender cambered wings has been reported in references 4 and 5. 

It should be noted that in the existing suction analogy method, as 
exemplified by reference 3, the edge suction forces predicted for attached 
flow are rotated so that they are normal to the cambered wing surface along 
the leading and side edges to produce the vortex lift effect. This is 
a direct extension of Polhamus' suction analogy originally developed for 
a flat wing (reference 6). However, experimental evidence (references 7 
and 8) indicates that the leading-edge vortex on a slender wing tends to 
migrate inboard as the angle of attack is Increased. This Implies that its 
suction force orientation depends on the local camber and the angle of attack, 
and is not always normal to the camber surface at the leading edge, as it 
is assumed in the existing method of suction analogy. Tlierefore, the migrating 
behavior of leading-edge vortex can not be predicted without modification of 
the current concept of suction analogy. In addition, the exact relation 
between the predicted thrust forces and edge suction forces has not been 
derived for a cambered wing. 
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The main purpose of this report is to present an improved method of 
suction analogy for slender cambered wings in subsonic and supersonic flow. 
The aforementioned deficiencies of the current method will be resolved, and 
comparison of experimental results with various existing methods for a 
variety of configurations will be given. 
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2. LIST OF SYMBOLS 


A 

b 

c 

c 

"a 

"av 


D. 

X 


D. . 

11 


C 

m 

C 

P 

‘"p(PM) 

AC 


c . 

tip 

T’ 

1, J , k 


K 


V, le 


K 


v.se 


aspect ratio 

span 

chord 

reference chord 

total axial force coefficient 

axial force coefficient due to leading-edge vortex 
sectional induced drag coefficient 
= C_ - (C I of non-cambered wings ) 

lc^=o 

total near-field induced drag coefficient 

total far-field induced drag coefficient in attached flow 

sectional lift coefficient 
total lift coefficient 
lift-curve slope at small a 

sectional pitching moment coefficient about Y-axis 
total pitching moment coefficient about Y-axis based on c 
pressure coefficient 

pressure coefficient calculated by Prandtl-Meyer theory 

lifting pressure coefficient 

sectional leading-edge suction coefficient 

sectional leading-edge thrust coefficient 

tip chord length 

unit vectors along X-, Y-, and Z-axes, respectively 
a unit vector normal to the wing leading edge (fig. 3) 
planform lift curve slope per radian at a = 0° 

leading-edge suction coefficient at one radian angle of attack 
side-edge suction coefficient at one radian angle of attack 
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Z-t- 3 + 


K 


a 

M 


N 

c 


an a-correction factor for the supersonic flow (eq. 54) 

Mach number, or number of integration points 

a unit vector normal to the wing surface 

a normal vector 

number of chordwise panels 



N 

s 


a unit normal vector to the wing surface at the leading edge 
number of spanwise strips on right wing 


n 

oo 


r 


s 



-V 

t 

CO 


U, V, w 

V 

V 

z 



X, y, z 


a unit vector normal to the freestream velocity vector 

streamwise distance of suction force center from the leading 
edge as defined in figure 7. Radial distance in figure 6. 

position vector 

local semi-span 

leading-edge suction force vector 
leading-edge thrust force vector 
a unit vector along the leading edge 

a unit vector along the freestream velocity vector 

induced velocity components along X- , Y- , and Z-axes, respectively 
velocity magnitude 

z- component of velocity in the vortex flow field (figure 7) 

circumferential velocity component in the vortex flow field 
(figure 7) 

rectangular coordinates with positive X-axis along axis of 
symmetry pointing downstream, positive Y-axis pointing to right, 
and positive Z-axis pointing upward 


<i> y^, 


a rectangular coordinate system obtained by rotating the XYZ 
system through an angle (f> about X-axis 


1 ’ 


'1’ 


^ 2 ’ ^ 2 ’ 


a rectangular coordinate system obtained by translating the 
X^Y^Z^ system along Y^-axis 

a rectangular coordinate system obtained by rotating the X^Yj^Z^^ 
system through an angle about Y^-axis 
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x,(y) 

z^(x,y) 

z^(x,y) 

a 

a 

Aa 

6 

A, A^, A2, 


x-coordinate of the leading edge 

camber surface ordinate measured from X-Y plane 

camber surface ordinate measured from the mean chord plane 

z-coordinate of the leading edge measured from X-Y plane 

angle of attack 

average local angle of attack including twist and camber 
angle of attack correction in supersonic flow (eq. 53) 
wing twist angle at y 
= sin 

M 

percent of elemental panel chord by which a control point 
on a leading-edge panel is moved downstream. See equations 
(47) and (48) 


6 

6 

c 

Y 

y. 

A 

A 


dx 

6 


deflection angle 
-1 dz 

= tan , c. 

ratio of specific heats (=1.4 for air) 
streamwise vortex density 
taper ratio 

leading-edge sweep angle 
dihedral angle 


= tan (Bz^/ay^) 
control surface 


elemental area vector 
-1 8z 

= tan , c. 

4x ^ 


Subscripts 


f 

.R.e, £ 

P 

r 


flap 

leading edge 

potential flow 
root 


se 


side edge or tip chord 



te 

vJle 

vse 

00 


trailing edge 
leading-edge vortex flow 
side-edge vortex flow 
freestream 



3 . THEORETICAL DEVELOPMENT 


For the present purpose, the wing is assumed to be thin, and is cambered 
and twisted. The flow field satisfies the linear compressible governing 
equation which is solved through the use of pressure doublets (references 1 
and 2) . To calculate the pressure distribution and other aerodynamic char- 
acteristics, the boundary condition of flow tangency must be satisfied. 

This condition and others will be developed in the following. 

3.1 Boundary Condition 

Assume that the wing surface (fig. 1) can be described by 
z = z^(x,y) (1) 

where * (x,y) is the ordinate of camber surface measured from the X-Y plane. 
Introduce a function f such that 

f = z - z^(x,y) (2) 

Therefore, a unit normal vector on the wing surface can be defined as 



The boundary condition on the wing surface requires that the total 
velocity component normal to the surface should vanish. Hence, 


[(V cosa + u)i + v^ + 


Using equation (3), equation (4) can be expanded to give 


[-V cosa 

CO 


3z 
c 

3x 



(5) 


where u is assumed small in comparison with cos a. 

To simplify equation (5), it is assumed that the pressure doublets are 
distributed on a mean chord plane which is defined to be a non-twisted plane 
inclined with a dihedral angle (()>) to the X-Y plane. Since conventionally 
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the airfoil camber is defined with respect to its chord line (denoted here 


as z^(x 2 , Y 2 ) and shown on fig. 2 ) and the twist angle is measured 

relative to the mean chord plane, it is convenient to express the camber 

slopes, and in terms of 4^, and (f) before equation (5) 

3x 3y “^2 “^2 ™ 

can be further simplified. This can be done through the following 
coordinate transformations. The original XYZ system is rotated about 
X-axis through an angle (j> (dihedral angle) , resulting in the 
system, and then is rotated through an angle (twist angle) about 

the axis to result in the system (fig. 1). In figure 1, the 

system is obtained from the X^Y^Z^^ system by a translation along 
the Y^ axis, just as the ^ 2 ^ 2^2 system is related to the ^ 2 ^ 2^2 

According to vector analysis, the results of such coordinate ro- 
tations can be obtained by a series of orthogonal transformations (p. 413, 

3z, 


ref. 9). It is shown in appendix A that 

dz 

-sin ^ cos OL 

tw dx2 

dz 

cos (b (cos a + sin ) 

^ tw dx 2 tw 

-sxn a + cos ot. , 

8 z tw dx„ tw 

c 2 


3x 


can be written as: 


9zc 

3x 


or, 


9x 


32 


( 6 ) 


(7) 


dz 

cos a + 3 — sin 

tw dx^ tw 


Similarly, -r — is given by 
dy 


3z 

c 

3y 


dz 


sin (b cos + dy„ 

^ tw 2 


cos 


dz 

c 

cos (b cos “ 3 — sin 

^ tw dy^ 


( 8 ) 


To simply equation (5) , note that the perturbed sidewash (v) is of the 


first order in magnitude. Following the thin wing theory, only first 
order terms will be retained in equation (5). Therefore, only the zero 
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3.2 Relation between Leading-Edge Thrust and Suction 

After the lifting pressure distribution is calculated, the leading- 
edge thrust coefficient Cj_(in the negative x-direction) can be determined 
by using the pressure distribution as described in references (1) and (2) . 
To calculate the suction coefficient (c^) from c^., the following steps 
are taken. 


Let ig be the unit vector along which the leading-edge suction force 
(1S|) acts in attached potential flow, t^ is a unit vector along the leading 
edge and is the unit vector normal to the wing surface at the leading 
edge. They are indicated in figure 3. It follows that: 




where n is given by equation (3) evaluated at the leading edge. To find 


t^, note that along the leading edge, equation (1) can be written as 
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Zj^(y) = z^(x,y) , X = 

The position vector of any point along the leading edge can be written as 
r(y) = x^(y)i + yj + Zjj^(y)k (13) 

The tangent vector is determined from 


dx. dz 

dr Z -t . ^ ? 

dy dy dy 


(14) 


where 


dy 

dy 


= tan A 

dz 3z 

= -r-^ tan A + ^ 
3 X 3y 


(15) 


(16) 


It follows that the unit tangent vector t^ is given by 

^ , I dr I 
dy ^ dy' 


(17) 


where 


= ‘ 

i 

/ 1 + tan^A + 

‘^^12 
My '' 


Substituting 

equations (3) 

and (17) 

into 1 




•t 

3 

-V 

k 

-?■ ">■ 

1 = n„ 

s i 

j-XJ 

az 

- c 

az 

- c 




3x 

3y 






dZj 



tanA 

1 

dy 


(18) 


/ 3z , 3z 9 

^ 'afl 


dz dz„ 3z dz. ^ 3z 3z 

= { ^ - 1] + j [ ^ ^ + tanA] + k[- — + tanA] } 

^ 8y dy 3x dy Bx By 


B z o 3 2 ry 

1 + ( — + ( — ^)^ !— I 

'■3x ^ '■3y ^ 'dy' 


(19) 


Since the suction force S is in the direction of i^ in attached potential 
flow and the thrust T is the component of S in the X-Z plane, it follows 
that 
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9z dz. 9z 


9z / / ^ « 9z „ “^ 

■.^tanA]>//l.(^)2.(^)2 |g 


'9y " 'dy' 
( 20 ) 


Thus, from equation (20), the sectional leading-edge suction coefficient 


(c^) can be obtained from the sectional leading-edge thrust coefficient 
(c^) as follows: 


c 

s 


/ o o 


I— I 

'dy' 


3 d z 2 

f + 1] ^ t- a 

d X d y 9 X 


92 3z . 

c , c ^ , ^2-,l/2 

+ ^ tanA] } 


( 21 ) 


where all quantities are evaluated at the leading edge, 9z^/9x, 9z /9y are 

• • d z 

given in equations (6) and (8), and £ is given in equation (16). 

dy 

3.3 Orientation of Rotated Suction Force 

In the method of suction analogy (reference 6) for a plane wing, the 
suction force predicted in the attached potential flow theory is rotated 
by 90° so that it is normal to the wing surface to simulate the vortex lift 
effect. Tf the same concept is used for a cambered wing, such as a delta 
wing with conically cambered leading edge (reference 7) , then the rotated 
suction force, being normal to the wing surface at the leading edge, would 
produce an increasingly large thrust component on the mean chord plane as 
angles of attack are increased. However, experimental data for a conically 
cambered wing (reference 7), reproduced in figure 4, indicate that at high 
angles of attack the locations of minimum pressure values will move inboard 
onto the planar portion of the wing. Since the experimental leading-edge vor- 
tex-induced suction force can be obtained by summing all vortex-induced suc- 
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tion pressure forces, the center of this suction force will also move inboard. 
Similar situation occurs on a delta wing with leading-edge vortex flap (ref. 8). 
This is seen in figure 5 from an examination of the axial force coefficient 



(C^) measured relative to the flat portion of the wing. The data show that 
stops increasing at about 16-17 degrees of angle of attack, indicating 
the leading-edge vortex has migrated over to the flat wing portion (refer- 
ence 8). Therefore, successful estimation of this suction force center is 
necessary for determining the appropriate orientation of the vortex force, 
as well as predicting the migrating behavior of the vortex. 

To be able to determine the location of this suction force center, a 
concept other than the suction analogy is needed. One possible way is to 
use the linear momentum principle of fluid mechanics. To get some idea on 
the flow field surrounding the leading-edge vortex, experimental results of 
Earnshaw (reference 10) are reproduced in figure 6. The velocity distribution 
was measured on a plane perpendicular to the freestream at a typical longi- 
tudinal station. Examining the circumferential velocity on a line a-a 
(fig. 6), its magnitude is approximately equal to the freestream velocity 
in a region from the vortex center to the leading edge, quadrant I, while 
it is less (about half of the freestream velocity in an average sense) 
inboard of the vortex center, quadrant III. The reason for the higher 
circumferential velocity outboard of the vortex center is probably due to 
the stronger effects of the vortex sheet and the "potential" flow. By 
"potential" flow, it is meant to be the non-separated flow component asso- 
ciated with the boundary condition being satisfied across the lifting surface. 
The results are quite similar in three other planes at 30, 50, and 70 percent 
of the root chord. It should be noted that the velocity shown in figure 6 
is the total circumferential velocity including both the potential and 
vortex flow effects. Since the circumferential velocity on line a-a inboard 
of the vortex center, quadrant III, can be assumed to be primarily due to 
the vortex effect, with the "potential" component being small, it is reason- 
able to assume that the vortex-induced normal velocity on line a-a is about 
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0.5 V . 

CO 

From the aforementioned observation and assumption, a proposed flow 

model for the present purpose is depicted in figure (7a), where is taken 

1 2 

as 0.5 V . Since c c ^ PV , by the suction analogy, is the vortex force 

03 g ^ 00 

per unit span, a width of unity for a section A-A is taken as shown in 
figure (7b) . Now, a control surface (a) is taken through the vortex center 
as shown in figure (7c). It is then assumed that the vortex force is equal 
to the vertical component of the force due to the momentum transfer through 
the control surface a. It follows that 

- CgC = J PV^(V • dT) (22) 

a 

From the assumption that the volume flow ratejv • dTjwill be: 

I 1 1 

|V ■ dx|= -V^ • 1 • r (23) 

Therefore , 

-c c^ P = I (PV - PV )V • dT 
s 2 oo J z . z 

a in out 

= -2( p i V • i V -1-r) 

or, 

r = c c (24) 

s 

Equation (24) is used to locate the assumed center of the vortex force 
and hence, the associated local slope for the orientation of the suction 
force. If the estimated suction force center is downstream of trailing 
edge, the corresponding camber slope for the force orientation is assumed 
to be zero. 

Although the present results for the location of vortex force center 
can be used to estimate the vortex contribution to pitching moment, unfor- 
tunately it is not possible to determine how the "potential" lifting pressure 
should be redistributed in order to estimate the "potential" contribution 
to pitching moment. Therefore, for the purpose of calculating pitching 
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moment, the vortex force will tentatively still be assumed to be at the 
leading edge, except that its orientation will be altered as discussed 
previously. 


3.4 Pressure Distribution and Total Forces and Pitching Moment 

The calculated longitudinal induced velocity (u) is assumed to be 
along the camber surface. Since in the linearized theory of both compres- 
sible and incompressible flow, the pressure coefficient (Cp) is related to 

u in the same expression: 
u 


-2 


C25) 


P V 

^ CO 

in the following the incompressible Bernoulli equation will be used to find 

a correction factor for high a in C . 

P 

The magnitude of undisturbed flow velocity on the mean chord plane is 
V^. By resolving this velocity along and normal to the camber surface, the 
components V cos (a + a 


tw 


6 ) and V sin(a + a - 6 ) can be obtained, 
c “> tw c 


as indicated in figure 8 for a = 0. It follows that the pressure coef- 

tw 

ficient (Cp) at any angles of attack can be obtained from: 

C' = 1 - [(V cos (a - 6 ) + u)^ + (V sin (a - 6 ) + w)^] 

P oo c oo c 


= -2 cos(a - 6 ) u/V = C cos(a - 6 ) 

C oo p c 

where 6 = tan ^(dz /dx„) and a = a + a .It follows that 

c ' c 2^ tw 

AC' = cos (a - 6 ) AC 
P c p 


(26) 


(27) 


where ACp is predicted by satisfying the boundary condition equation (10a) . 
Physically, equation (27) Implies that the predicted vortex strength from 
the attached potential flow theory is allowed to interact with the free 
stream velocity component tangent to the camber surface to produce the lifting 
pressure (AC^) (reference 6) . 
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For a cambered wing with or without dihedral, additional lifting 
pressure will be generated from the interaction of a freestream component 
with the strearawise vortex density (y^) > shown in figure 9. Adding this 
component of lifting pressure to that given in equation (27) results in: 


AC" = cos(a - & )AC - 2y sina sin($ + d)),a = a+a^ 
p _ c P X y tw 

-1 


(28) 


(29) 


where may be defined as "local dihedral" angle and may be different 
everywhere on the cambered wing. On the other hand, the geometric dihedral 
angle ((> is constant along y for a given cambered or noncambered wing panel. 
For a noncambered wing, = 0 but (|) may not be zero. Note that in the 
normal lifting condition, is positive; because when viewed as a vector, 
it is pointing in the downstream (i.e., positive X-axis) direction. 

The lifting pressure (AC^) is acting normal to the local camber surface, 
i.e., in the n direction. This pressure force will be decomposed into n^ 
and directions to determine the lift and drag components, respectively, 
where (see fig. 10) 


sin a i + cos a k 


cosai+sinak 


(30) 

(31) 


Using equation (3) for n and then determining n • n^ and n • t^, it can 
be shown that 




'd,p 


_ 1 


dz 


(ax' 

c 

J ^^P 


^He 




dz 

_ 1 
C 

/ 

C 

“ ax 


ie 



sina + cosa) 


3z „ 3z „ 

^ ( 32 ) 


3z 3z 

3 3- (^) + dx (33) 


To find the pitching moment, the j-component of 
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3z 3z 

- c -t c -t ->- 

, A 3x ^ 3y ^ 

(Axi + Azk) X / ^ 9 ^ o ' 

/l + + (^)2 

is needed, where Ax and Az are the moment arms of the pressure elements 
referenced to the coordinate system. It follows that 

dz 


m,p ^ 


/ 


te 


AC 


& 


3z 


(Ax + Az — ) jjl + (— ) + (— ) dx (34) 


"Jle 


The contribution of leading-edge vortex effect to the sectional aero- 
dynamic characteristics follows the same types of expressions given in 
equations (32)-(34). Hence, 

3z 


~l,vZe 


'-s 


3z 

'^d,vAe " ^ " IT 


- 0 “ “ (Ax„ + Az, 

m,vZe c Ze i 

3z 3z 



(35) 

(36) 

(37) 


c c 

where -r — and are to be evaluated in accordance with the discussions 

3x 3y 

of Section 3.3, regarding vortex movement. 

Finally, spanwise integration of sectional characteristics given in 
equations (32)-(34) and equations (35)-(37) will produce the total force 
and moment coefficients: 


2 r b/2 

'L,P 


c dy 

(38) 





2 f 

'D,p “ 

o 

‘^d,p 

c dy 

(39) 

= T. 

’m , p 

c 

m,p 

2 , 
c dy 

(40) 

0 rb/2 



'L.vJie SJ 
o 

^z 

„ dy 

jVJle 

(41) 
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2 1 

..b/2 

'DjVAe S j 

Cj „ c dy 
d jVJle 
o 

. _ 2 

.b/2 

'm,vJle Sc 

j "m,v£e " 

o 


( 42 ) 


(43) 


The above integration can be performed by applying the trapezoidal rule, 
as shown in reference 11. 

The side-edge contribution to the vortex-lift effect can be obtained 
by resolving the side-edge suction force, which is normal to the camber 
surface along the side edge, into n^ and t^ directions to determine lift 
and drag components, respectively. They have similar forms as those 
in equations (35)-(36) . The pitching moment contribution follows 
equation (37) closely. Thus, 

dz 


- 


LjVse 


„ (-;: — sin a + cos a) 

se '9x 


- ‘^tip 


D ,vse 


m,vse 


^ f c 

J S , 

r dz 

^ r c (- -^ 

j '-s,se 3 x 


o 


cos a + sin 
& 

(Ax + Az ^) 

se se se 3x 


f 


dx (45) 


/ ^9 8z „ 

^ + (- 3 T) 


where the integration is performed over the tip chord and c^ is the local 
side-edge suction coefficient. The determination of side-edge suction 

3z 

coefficient c is discussed in references 2 and 12. The values of r — 
s,se 3x 

3z 

and ^ — in equations (44)- (46) are determined along the tip chord without 
oy 

assuming inboard movement of side-edge vortex. This is based on the 

dz 

assumption that on usual configurations does not vary too much near the 
tip and c^ is usually much smaller than c^ along the leading edge so 
that the vortex center will be quite close to the side edge. 

In appendix B, analytical expressions to describe the geometry of 
leading-edge plane flaps and conical camber are derived. A numerical 
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method for calculating 
in appendix C. 


3z 

gy from the prescribed camber ordinates is given 
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4. NUMERICAL RESULTS AND DISCUSSIONS 

The computer program used In references (1) and (2) serves as a basis for 
further development to account for all features developed in Section 3. The 
three-dimensional effects (called sweep effect in references (1) and (2)) on 
choice of control-point locations have been re-examined. Since large geometric 
dihedral angles are now allowed in the program, the accuracy of calculation for 
configurations with large dihedral in both attached flow and separated flow 
must be assessed. The results of this re-examination and assessment, along 
with cambered wing data comparisons will be presented and discussed next. 

4.1 Three-Dimensional Effects on Control-Point Locations 

In the present constant pressure panel method, the control point on 
an elemental panel is chosen based on a two-dimensional consideration to 
match pressure values with those given by the exact thin airfoil theory. 

However, it is well known that for a wing with large sweep angles and/or 
small taper ratios, the outboard pressure distribution can be quite different 
from a two-dimensional shape. To account for this effect, it was proposed 
in reference 2 to move downstream the control points of the leading-edge 


panels by an 

additional amount in percent. A, of the 

elemental panel chord 

where. 





for M < 

1, 

> 

t| 

O 

for A < 60^ 

(47a) 



= 0.35(A-60) 

for A > 60° 

(47b) 

for M > 

1, 

A, = 3.0 

for A < 60° 

(47c) 



1 





= 0.35(A-60) 

for A > 60° 

(47d) 


A = 

A. + A- 


(47e) 
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Equation (47e) appeared to be accurate for highly swept wings (A > 70°) 

with any taper ratio examined in reference 2. However, for moderately swept 

wings (A < 70°) with zero taper ratio, the results may not be accurate. This 

is illustrated in table 1(a). The assessment of accuracy is based on the 

ratio of near-field to far-field induced drag coefficients, /C^ , in the 

i ii 

attached flow solution (i.e. with 100% leading-edge thrust) (ref erence 13). 

This ratio should be as close to unity as possible. A value of less than 1.0 

for Cjj /Cp implies that the predicted leading-edge suction is too large, 
i il 

This appears to be true for all cases in table 1(a) , except case (2) . 

Apparently, not only the taper ratio (see cases (1) and (3)), but also the 
aspect ratio (see cases (1) and (5)) will affect the accuracy of the attached 
flow solution. To improve the accuracy, the following empirical factor, 
in percent, was used for all M's: 

A 3 = 3.5 (1-A)^ ^ (48a) 

A = A 3 + greater of (A 2 , A 3 ) (48b) 

The aspect-ratio factor in equation (48a) was chosen such that for A = 2.0, 
the aspect-ratio factor is one. The results are presented in table 1(b) and 
appear to have significant improvement. Therefore, equation (48b) will be 
used in all following results. It should be noted that the adjustment 
represented by equation (48b) is applied to the outboard region of the wing 
only. For inboard regions, equation (47e) is still applicable. 

4.2 Effects of Dihedral Angles on Flat Wings 

Effects of large dihedral angles on wing aerodynamic characteristics 

have been experimentally investigated in reference 14 with V-tails. In 

figures 11(a) and 11(b), the predicted results for two non-cambered V-tails 

in attached flow are compared with data. These tails are called in reference 14 

Tail A and B and have aspect ratios of 5.55 and 3.70, respectively. Results 

by the quasi-vortex-lattice method (QVLM) (ref. 11) are also shown for the 

largest dihedral angles presented. The agreement between the present results 
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Table 1. Predicted Lift Factors and Induced Drag Ratio 
for Various Flat Configurations 


(a) Present Method with equation (47e) 



Geometry 

M 

K 

P 

Sl,ie 

^,se 

S/S.. 
1 11 

Methods 

N xN 
c s 


Delta 


3.16613 

3.88620 

0. 

0.684 

Present 

10 

X 

12 

(1) 

A=50°, X=0. 

0. 

3.03177 

3.18721 

0. 

1.087 

ref. 11 

7 

X 

12 


Delta 


1.54912 

2.95563 

0. 

1.0958 

Present 

10 

X 

12 

(2) 

A=74°, A=0. 

0. 

1.43638 

2.93962 

0. 

1.0815 

ref. 11 

7 

X 

12 


Delta 


2.94631 

3.31705 

0.20858 

0.803 

Present 

10 

X 

12 

(3) 

A=50°, A=0.1 

0. 

2.86825 

2.90379 

0.20250 

1.0423 

ref. 11 

7 

X 

12 


Delta 


2.54488 

3.89359 

0. 

0.7755 

Present 

10 

X 

14 

(4) 

A=63.43°, A=0. 

0.7 

2.39327 

3.12126 

0. 

1.0814 

ref. 11 

10 

X 

14 


Delta 


4.90702 

4.49280 

0. 

0.8661 

Present 

7 

X 

15 

(5) 

A=20°, A=0. 

0. 

4.84672 

4.28371 

0. 

1.0615 

ref. 11 

7 

X 

15 


A= 45° 


1.63405 

1.74318 

0.38427 

0.9216 

Present 

7 

X 

15 

(6) 

A=2. 

0, A=0.5 

0. 

1. 62226 

1.68169 

0.37001 

1.0087 

ref. 11 

7 

X 

15 



^0° 


1.83039 

3.57093 

0. 

0.8179 

Present 

9 

X 

14 

(7) 


_^65° 

0. 










same L.E. 
length in 
both regions 


1.66340 

3.03010 

0. 

1.1199 

ref. 11 

9 

X 

14 


* 

Attached flow results. 
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Table 1. (Concluded) 


(b) Present Method with equation (48b) 



Geometry 

M 

K 

P 


^,se 

1 11 

Methods 

N xN 
c s 



Delta 


3.15358 

3.19500 

0. 

1.135 

Present 

10 

X 

12 

(1) 

A=50°, X=0. 

0. 

3.03177 

3.18721 

0. 

1.087 

ref. 11 

7 

X 

12 


Delta 


1.54912 

2.95563 

0. 

1.096 

Present 

10 

X 

12 

(2) 

A=74°, X=0. 

0. 

1.43638 

2.93962 

0. 

1.0815 

ref. 

11 

7 

X 

12 


Delta 


2.93790 

2.82564 

0.20593 

1.113 

Present 

10 

X 

12 

(3) 

A=50°, X=0.1 

0. 

2.86825 

2.90379 

0.2025 

1.0423 

ref. 

11 

7 

X 

12 


Delta 


2.53932 

3.24207 

0. 

1.056 

Present 

10 

X 

14 

(4) 

A=63.43°, X=0, 

0.7 

2.39327 

3.12126 

0. 

1.0814 

ref. 

11 

10 

X 

14 

(5) 

Delta 

0. 

4.88981 

4.14896 

0. 

1.261 

Present 

7 

X 

15 

A=20°, X=0. 

4.84672 

4.28371 

0. 

1.0615 

ref . 

11 

7 

X 

15 

(6) 

A=45° 

0. 

1.63218 

1.68857 

0.38310 

1.008 

Present 

7 

X 

15 

A=2.0, X=0.5 

1.62226 

1.68169 

0.37001 

1.0087 

ref. 

11 

7 

X 

15 


^80° 


1.82592 

3.10805 

0. 

1.055 

Present 

9 

X 

14 

(7) 

^65° 

same L.E. 
length in 
both regions 

0. 

1.66340 

3.03010 

0. 

1.1199 

ref. 

11 

9 

X 

14 


5? 

Attached flow results. 
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and those predicted by the QVLM Is very good. The experimental data are 
all reasonably well predicted. 

A nonplanar configuration exhibiting edge-separated vortex flow is 
shown in figure 12. The predicted aerodynamic characteristics are compared 
with experimental data (ref. 15) in figure 13. Several points are of 
interest . 

(1) If the vortex lift is assumed to exist along all edges (i.e., full 
vortex lift curves), overall lift and drag are over-predicted. This has 
been verified in reference 15 with the conventional vortex-lattice method. 

(2) If the outboard panel (identified as section II in fig. 13) is 
assumed not to carry vortex lift, but the leading-edge suction is also lost, 
the predicted results in and agree well with experimental data. This 
assumption is consistent with experimental observation that the organized 


vortex flow is weak on the outboard panel. 

(3) The trend of variation due to tip dihedral is correctly predicted, i.e., 

increasing the tip dihedral reduces the C, and C,.^, and makes C more positive. 

L u tn 

dz 

4.3 Effect of on Leading-Edge Suction Coefficient for Cambered Wings 
Among the cambered wings examined so far, the delta wing with conical 
camber tested in reference 7 appears to be the most critical case in that 


along the leading edge the camber slopes amount to 0.3845 and -1.4402 for 

and ' ^^c , respectively. The camber geometry used in the present calcula- 

dx2 dy^ 

tion is derived in appendix B. From numerical calculation, it was found that 


a large negative value of c can greatly reduce the leading-edge suction 

9z 9z 

coefficient when equation (21) is directly used. (Note that -r — , -r — and 

dz dz 

- 3 — in eq. ( 21 ) are related to c and c in eqs. ( 6 ), ( 8 ), and (16), res- 


pectively) . This is illustrated in figure 14. It is seen that without 
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dz 

(> 

setting - 1 — to zero in equation (21), c c is too low and hence the lift is 
uy 2 s 

underpredicted. This is perhaps because the boundary condition (eq. 10b) 

is applied on the mean chord plane and the doublet is distributed also on 

the same plane, instead of the actual camber surface. Because of the large 
dz 

magnitude of c near the leading edge for this conically cambered wing, 

neglecting it will violate the small perturbation assumption in a small 

region near the leading edge. However, the present method is consistent 

d z 

with the usual practice in the thin wing theory. This means that c has 

dyz 

not been used in calculating the doublet strength. Therefore, in all calcula- 

d z 

tions to be presented below, c will be set to zero in equation (21) . 

dy2 

4.4 Effect of Vortex Force Orientation 

Again, the conically cambered slender wing of reference 7 will be used 
to illustrate the effect of vortex force orientation. By setting r=0 in 
equation (24) , the vortex force will then be normal to the camber surface 
along the leading edge, an assumption being used in the existing methods. 

Figure 15 Indicates that the present proposed theory (r = c c) can signifi- 
cantly improve the predicted results. Therefore, all following calculations 

will correspond to r = c c. 

s 

As Indicated in Section 3.3, the relation r = c c is mainly used to 

s 

calculate the location for the vortex force orientation on a cambered wing. 

It is of interest to compare the location of vortex force centroid predicted 
by the present method with the leading-edge vortex location given by other 
methods. Figure 16 presents a comparison with the results by Brown and Michael 
in reference 16. The present method predicts the vortex force centroids 
slightly outboard of the vortex location given by reference 16. The model 
of Brown and Michael has further been employed in reference 17 where theoretical 
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results were compared with data for different planforms. The comparison 
indicated that the vortex force centroid shown in figure 16 by the present 
method and the vortex location by ref. 16 is too much outboard of the ex- 
perimental vortex location. The same conclusion is true when compared with 
Smith's calculation (ref. 18). It should be noted that one additional 
effect included in reference 18, but not modelled in references 16 and 17, 
is the leading-edge free vortex sheet which is to increase the loading near 
the leading edge. This Implies that the vortex force centroid will not 
coincide with the vortex location, but will be outboard of it. Therefore, 
the present method should be used only to find the vortex force centroid, 
but not the concentrated vortex-core location. 

4.5 Cambered Wings in Subsonic Flow 

Several cambered slender wings have been used to correlate the pre- 
dicted results with data. The results for the aforementioned conically 
cambered delta wing are compared with those by the vortex lattice method, 
VLM-SA, (ref. 3) in figure 17. The present results appear to agree better 

with data than those by VLM-SA in AC and C . However, the difference in 

D Lj 

the predicted C^ is smaller than expected, despite of the fact that in 

VLM-SA the vortex force is assumed to be normal to the camber surface at 

the leading edge (see figure 15 for the effect) . One possible reason is 

3 z 

that in reference 3, c is taken to be zero in the denominator of equations 

3y 

(32) through (37), resulting in Increased C and AC . The predicted AC 

Lj iJ JJ 

at high angles of attack by the present method are slightly too low. Whether 

this is due to additional separation at the juncture of conical camber and 

the plane wing section is not known. On the other hand, the present pitching 

moment is too negative at high lift coefficient. 

A highly cambered wing was tested by Squire and reported in reference 19. 

This wing has camber slopes, ^^c and ^^c , along the leading edge as high as 

3x 3y 



0.45 and -2.1, respectively. From the comparison made in figure 18, , 

and are all well predicted by both the present method and VLM-SA. 

The jump in experimental at a=20° was explained in reference 19 to be 
due to merging of two vortices - one being the leading edge vortex and 
the other one being inboard vortex due to camber shoulder separation. 

Recently, the use of leading-edge vortex flaps to increase the lift- 
drag ratio under maneuvering conditions has been of great interest. Some 
low speed results for a 74°-delta wing have been reported in reference 8. 

The geometry of the configuration used in the present method is shown in 
figure 19, where the dashed lines form the altered boundaries to define the 
planform in the present and VLM-SA computer programs. The undeflected 
trailing-edge flap on the test model is ignored for simplicity. The results 
are presented in figure 20. Note that all coefficients are based on total 
planform area excluding the flap. It is seen that the predicted results 
agree well with experimental data, as do those for the VLM-SA method. Note 
that the ability of the present method to predict the vortex force orienta- 
tion makes it possible to estimate when the leading-edge vortex will move 
onto the plane wing. As shown in figure 21, the magnitude of vortex axial 
force coefficient starts to drop at a=16° , implying that for a > 16° 

the vortex will migrate onto the wing. In reference 8, the a at which this 
migration occurred was estimated to be 17°. The present results are there- 
fore seen to be in good agreement with this reference. 

Another configuration with full-span leading-edge flap was tested in 
1958 at subsonic and supersonic speeds (ref. 20). The wing-body combination 
consists of a delta wing of aspect ratio 2.0 and NACA 0003-63 airfoil in the 
streamwise direction. The test results at M=0.7 are compared in figure 22 
with the predicted values for a configuration with the leading-edge flap 
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deflected down by 15 degrees normal to the hinge line. It is seen that the 

theoretical results by the present theory and VLM-SA agree quite well with 

the data. It should be noted that the fuselage in reference 20 is not 

modelled in either the present theory or the VLM-SA. 

Two highly cambered wings designed by the method of vortex- lattice 

and suction analogy have been tested in subsonic flow. The results are 

presented in figures 23 and 24 together with theoretical prediction. It is 

seen from figure 23 that the present method offers significant improvement 

over VLM-SA in drag prediction. Presumably, this is due to the more accurate 

vortex force orientation calculated in the present method. For the cranked 

cambered wing shown in figure 24, the present predicted drag level is slightly 

too high at low C . The reason is not known. For C and C , the present 
i-i Lt in 

results are in reasonably good agreement with data, except that the for 
the cranked cambered wing is underpredicted at moderate angles of attack. 

These results are not, in general, any better than those of the VLM-SA with 
the exception of the prediction. 

4.6 Wings in Supersonic Flow 
Non-Cambered Wings 

The method of suction analogy was first applied to delta wings in super- 
sonic flow by Polhamus (ref. 22). If the concept is directly used on a delta 
wing of A = 4/3, results shown in figures 25 and 26 can be obtained. Note 
that the dashed curves denoted as the present theory without correction are 
the results by direct application of Polhamus' concept. So are the curves 
marked as original suction analogy. It is seen that the theory overpredicts 
C. and underpredicts AC slightly. Examination of the test cases in reference 

1j U 

22 indicates that the same trend persists on a delta wing of A = 1.0 presented 
in that reference. Although the attached flow theory (with 100% leading-edge 
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suction) tends to predict C at low a rather well, the drag is underpre- 

Xj 

dieted. In addition, the test model used in reference 23 has sharp edges 
so that the concept of partial leading-edge suction Investigated in 
reference 24 is not applicable. In an effort to improve the situation, 
it is helpful to study the flow field on a 65®-delta wing reported in 
reference 25. It was indicated in reference 25 that in supersonic flow 
the leading-edge separation starts at a higher angle of attack than that 
in subsonic flow on the same wing. The vortex flow region is diminished 
quickly as the Mach number is increased at a given angle of attack. In 
addition, the separated flow characteristics on a wing with small round 
leading edge (RAE 101 section) are quite similar to those with sharp 
leading edge (biconvex section). In fact, the longitudinal aerodynamic 
characteristics with fixed transition are practically the same for both 
airfoil sections. Since there is no doubt about the existence of separated 
vortex flow on slender wings in supersonic flow, the discrepancy between 
the results by current suction analogy and data must be nonlinear effect 
not accounted for in the linear potential flow theory. Of course, this 
conclusion is based on the assumption that the current suction analogy is 
still applicable in supersonic flow. 

To find a correction to the linear potential flow theory, it is noted 

that the flow expansion on the suction side of the wing makes the vortex 

flow separation possible. Therefore, it is of interest to compare the 

suction pressure predicted by the linear theory with that given by the Prandtl- 

Meyer theory. In the two-dimensional case, it is shown in reference 26 (p. 386) 

that the linear theory predicts always a higher suction pressure. For a flat 

plate at an angle of attack a, the upper surface pressure coefficient (C ) 

P 

given by the linear theory is (ref. 26, p. 386) 



On the other hand, the Prandtl-Meyer theory shows that (ref. 26, p. 383) 




, sin B . y-1 
~ cos 26 '' 


(- 


- 2„ 

sin 3 


Y 

a-1 


- 1 


Y - cos 23 


(50) 


where 3^ = sin ^ ^)and 3 = sin ^ ^) which must be found from (ref. 26, p. 377) 


-a + 3 


- tan ^ tan 3) = 3 tan ^ (/^^ tan 3 ) (51) 

yj y~i V Y“i “ yy~i VY“1 “ 


To make = Cp(PM), a correction (a') to a in equation (49) must be made 
such that 


2 (g + g') 


J 


M - 1 


"p(PM) 


or. 




" ”2 %(PM) 


(52) 


Equation (52) must be further modified based on the following two consider- 
ations. Firstly , equation (52) has been obtained by considering the upper 

surface pressure coefficient only, not the lifting pressure (AC ). Assuming 

P 

that the lower surface compression does not affect the upper surface vortex 
flow so that the linear theory prediction does not require modification, 
then the corresponding angle of attack correction to AC^, Aa, will be given by 

4 (g + Aa) _ 2g 2 (g + g' ) 

- 1 - 1 /m^ - 1 

V 00 V 00 \/ CO 

or, 

A 1 . 

Aa = 

Secondly , experimental evidence (ref. 25) indicated that as the Mach number 


is increased at a given a, the vortex flow region will quickly diminish faster 
than predicted by reference 22. This means that a' should decrease 



in magnitude as is increased. Presumably this is due to the three- 
dimensional effect, effect of conical shock above the vortex flow (ref. 25), 
etc. One simple way to account for these effects is to assume that a' varies 


[2 

inversely withjM^ - 1, in much the same way as does (see eq . 49). A 

Mach number of 1.4 will be chosen as a reference (see below), because it 
is usually regarded as the lower bound in Mach number for the linear super- 
sonic theory to be applicable. Combining the above considerations, the 
a- correction, Aa, can now be written as 

Aa - (-C. (53) 

where 

M^<_ 1.4 (54a) 


K =i 




4-1 


M - 1 


M >1.4 


(54b) 


Equation (53) is illustrated in figure 27 with and without equation (54b) 
for a = 10 deg. The magnitude of Aa for different M^ and a by equations 
(53), (54a) and (54b) is indicated in figure 28. 

The results of applying Aa correction to the calculation for a delta 
wing of aspect ratio 4/3 have also been presented in figures 25 and 26 for 
M^ = 1.4 and 1.8. It is seen that a significant improvement in prediction 
is apparent. 

In reference 22, the results by the original suction analogy method 
were compared with data from reference 27 for a delta wing of aspect 
ratio of 1.0 at M^ = 1.97. This comparison is reproduced in figure 29. 

It is seen that the lift is overpredicted and the drag at a given lift 
is underpredicted by the original method. On the other hand, the present 
modified method predicts both the lift and the drag quite well. 
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A slightly cropped delta wing of 76° sweep angle has been tested 
in supersonic flow with results reported in reference 28. The plane 
wing had 4-percent-thick circular arc airfoils. The test data are compared 
with the theoretical prediction in figure 30. In the present theory, if 
the side-edge (SE) vortex lift effect is included, the lift will be slightly 
overpredicted (about 4% at a = 20 deg.) and the drag slightly underpredicted. 
The pitching moment is slightly more negative at high than the data show. 
On the other hand, if the side-edge vortex lift effect is ignored, good 
agreement between theory and data can be achieved. Whether this is 
true for wings with higher taper ratios is not known, since most available 
data were obtained on wings with zero or rather small taper ratios. The 
results by the Middleton-Carlson linear theory are also shown in figure 30. 

It underpredicts the lift at high a's and overpredicts the drag slightly. 

Finally, the results for the cropped plane delta wing of 65°- sweep 
angle reported in reference 25 are presented in figure 31. The fuselage 
is not modelled in the present program. It is seen that the theory with 
full vortex effect overpredicts by about 10% at a = 10 deg. Excluding 
the side-edge vortex lift effect does not improve too much the agreement. 

One possible reason for the discrepancy is the presence of a fuselage 
which has a diameter- to-wing span ratio of 0.2083, and has a forward 
portion of 0.75b ahead of wing root chord. Squire (ref. 25) observed 
that the fuselage tends to increase the angle of attack at which the 
vortex flow becomes significant. The main cause is the conical shock from 
the fore body. From figure 31, it is shown that the present theory with 
the potential flow component only (0% leading-edge suction) seems to 
predict the results rather well. The good agreement in this case may 
be fortuitous. 
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Further investigation of the fore body effect in supersonic flow is needed. 


Cambered Wings 

For cambered wings, equation (53) should be modified by using some 
average local a: 


Aa 


K (-a 
a 





(55) 


where is still given by equations (54). In the following calculation, 
a is assumed to be the arithmetic mean of local a including twist and 
camber at all control points. 

Several configurations of leading-edge flaps have been tested at = 1 
and 1.8 and reported in reference 23. Results for deflection angles of 
36.45 degrees and 60 degrees are presented in figures 32 and 33. It is 
seen that the lift is well predicted in all cases, in particular at = 1,8 
At a transonic Mach number of 1.4, the predicted is slightly higher, 
with discrepancy being 8% or less. However, in all cases the drag is 
underpredicted for positive C, and overpredicted for negative . The 
same situation has been demonstrated in reference 30 (fig. 4C of that 
reference) for a configuration with leading-edge flaps deflected 20 deg. 
Possible reasons for this are flap hinge separation at positive as 
indicated by the author in reference 23 and the nonlinear effect on wave 
drag. At negative lift, the drag data fall between the results for 
attached flow (100% leading-edge suction) and vortex flow. Whether 
this implies partial leading-edge vortex flow on the lower surface is 
no t known . 

To avoid the effect of possible leading-edge flap separation, a 
smoothly cambered wing would be more suitable to test the theory. Data of 
such a wing are available in reference 28 for a wing designed to cruise at 
= 3.5 and design of 0.1. The results are presented in figure 34. 
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Again, the present theory without the side-edge vortex lift effect appears 
to be more accurate for C and AC • However, it gives less accurate 

J-i jJ 

results for C at high . 

in Jj 
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CONCLUDING REMARKS 


An aerodynamic panel program has been developed for cambered wings 
exhibiting edge-separated vortex flow at subsonic and supersonic speeds. 
Semi-empirical but simple methods have been developed to determine the 
appropriate position and hence orientation of vortex force on cambered 
wings, as well as to improve the prediction in supersonic flow. Com- 
parison with available data indicates that (1) the program is accurate in 
predicting subsonic aerodynamic characteristics of highly cambered wings 
with vortex separation effect. The improvement in drag prediction over 
currently available methods can be significant. (2) The program is also 
accurate in predicting supersonic aerodynamic characteristics of slender 
wings with vortex separation. This has been achieved by employing an effec- 
tive angle of attack in the suction analogy method. 

In view of the improvement made in the present prediction method and 
possible further refinement, it would be of interest to: 

(1) use the program to design a supersonic wing so as to make 
efficient use of the vortex flow; 

(2) examine the fuselage effect on wing vortex flow more closely; 

(3) reexamine the drag associated with using leading-edge flaps in 
supersonic flow; 

and (4) make additional verification studies in which the significance 

of side-edge vortex lift in supersonic flow can be investigated. 
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APPENDIX A 


Derivation of Expressions for Camber Slopes, 


3z 

and 


c 




9z 9z 

As indicated in Section 3, the camber slopes, c and c, based 

3x 3y 

on the original XYZ system (see fig. 1) can be expressed in terms of 
camber slopes based on the local chord line. To achieve this, two rota- 
tions of the coordinate system are needed. The first is about the 
original X-axis through a dihedral angle (^) , and the second is about 

the new Y-axis (Y-,-axis) through a local twist angle (a ). The resulting 
J. tw 

expressions for ^^c and ^^c are derived below. 

3x 3y 

3z 

c 

A. 1 Calculation of — 

3x 

Consider a vector with components (dx, dy, dz^) tangent to the camber 
surface at some point, P. This vector has the component PQ in the X-Z plane 
(fig. 2). According to vector analysis, the vector (dx, dy, dz^) based 
on the XYZ system can be expressed in the X^Y^Z^ system (dx 2 , dy 2 , dz ) 
through the following relation (p. 413, reference 9) : 




COS a 

tw 

0 

-sin 

tw 


~i 

0 

0 



^dx \ 

dyz 

= 

0 

1 

0 


0 

cos (}> 

sin 

((> 


dy 



sin a 

tw 

0 

cos a 

tw 


0 

-sin (j) 

cos 



dz 

\ 7 


COS a 

tw 

sina 

tw 

siiKp 

-sina 

tw 

costj) 


/ dx ^ 

0 

COS (f) 


sin 

<!> 


dy 

sin 

tw 

- sin(j) 

cosa 

tw 

cosd) cosa 

tw _ 


{“‘J 


It follows from equation (A.l) that 


dx„ = dx cos + dy sin sin * - dz sin a. cos d> 
2 tw tw c tw ^ 

dy 2 = dy cos <() + dz^ sin (|) 

dz = dx sin a - dy sin tp cos a + dz cos d) cos a 
c tw ^ tw c ^ tw 


(A.l) 


(A. 2) 
(A. 3) 
(A. 4) 
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3z 


dz 


To calculate in terms of 


(A. 4). Hence, 


y=constant 


dx„ = dx cos a. - dz sin a cos ij) 
2 tw c tw 

dz = dx sin + dz cos (|) cos a 

c tw c tw 


, set dy = 0 in equations (A.2)- 

(A.5) 

(A. 6) 


Equations (A. 5) and (A, 6) can be solved for dx and dz^ through Cramer s rule 
to give: 


dz = 
c 


-dx„ sin + dz cos a 
2 tw c tw 


dx = dx„ cos a. + dz sin a 

2 tw c tw 


3z 


(A. 7) 
(A. 8) 


To find equation (A. 7) is divided through by equation (A. 8) to yield: 

3x 


^ -sin a + c cos a 

3z tw tw 

c _ dX 2 

3x — — 


(A. 9) 


y . dz • \ 

cos (f)(cos ct^ + c sin ) 

^ ^ tw -z tw 

^^2 


or. 


3z 
c 

3x 


, dz 

-sin + c cos a. 

tw tw 

dx2 

. dz 

cos a + c sin 

tw tw 

dX2 


(A. 10) 


3z 


dz 


For 0 = 0, let = tan 0 and 


tan 6 . Both 0 and 6 are 
c c 


indicated in figure 2. Equation (A. 10) then becomes 

cos5 + sin6 cosa 


tan 0 


tw 


tw 


cosa cos5 + sin6 sina^ 
tw c c tw 


sin (6 - ) 

' c tw 

cos (S - ) 

c tw 


= tan (6 - ) 

c tw 


Hence, 

0=6 “ u (A. 11) 

c tw 

The result in equation (A. 11) is seen to be correct from a geometric point 
of view, because in figure 2 , 0 is negative and 6^ is positive. 
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A. 2 Calculation of -r — 

di 

c c 

To calculate in terms of , set dx = 0 in equations (A. 2) -(A. 4). 

It follows that: 


dy 2 = dy cos <() + dz^ sin cj) 

dz = -dy sin d) cos + dz cos d> cos 
c tw c tw 


(A. 12) 
(A. 13) 


Solving equations (A. 12) and (A. 13) for dz^ and dy through Cramer's rule, 
it is obtained that 


dy« sin d) cos a + dz cos 


dz = 
c 


dy„ cos <t> cosa^ - dz sin 
-'2 ^ tw c 


3z sin d) cos + - 3 — cos 

c _ tw dy? 

3y ~ dz 

cos <p cos a - — sin 

tw dy - 


(A. 14) 
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APPENDIX B 


Geometry of Leading-Edge Plane Flaps and Conical Camber 
In this report, leading-edge plane flaps and conical camber have been 
used extensively to demonstrate the applications of the present program. 

These two kinds of aerodynamic devices have heen represented analytically 
in the program. For the convenience of readers, these analytical expressions 
are developed in the following. 


B.l. Leading-Edge Plane Flaps 

A leading-edge plane flap, as shown in figure B.l, is defined by 
four corner points, 1-2-3-4. In applications, points 1 and 2, or points 
3 and 4, may coincide. Therefore, to define a normal vector to the plane, 

It xs convenxent to use the two vectors and given by: 

h “ ^^4 “ ^1^^ ^^4 " ^1^^ ^^4 " ^1^^ 

^2 ^^3 “ ^2^^ ^^3 “ ^^3 “ (B.2) 

A normal vector (ij) to the plane flap can now be constructed by using T^ and T 

= T^ X T^ = 


j k 




^4 - 

X3 - 

^1 

X2 

1 

= 1 

[(Y 4 - Ti) 

(Z3 - 

Z2) 


[(X3 - X2) 

(24 - 

Zi) 

+ic 

[(X4 - x^) 

(Y 3 - 

T2) 


= Ni + Nt + N'S 

X y z 

It follows that if the equation 
ax + by + cz + d = 0 


T4 - yi ^4 - ^1 
^3 - ^2 ^3 ■ ^2 

- (T3 - V2> (^4 - 23 )] 

- (x^ - x^) (Z3 - Z2>] 

- (X 3 - X 2 > (y^ - y^^)] 


(B.3) 

of a plane containing the flap is written as 

(B.4) 
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the constants, a, b, and c can be chosen as follows: 


a = N, b = N, c = N 
X y z 

3z 


(B.5) 


3z 


c c 

The camber slopes, and , needed in determining the aerodynamic 

characteristics can then be obtained by differentiating equation (B.4) after 
z is replaced by z^: 


3z 
c 

3x 

3z 

c 

3y 


a 

c 

b 

c 


N 

_ 

N 

z 


(B.6) 


(B.7) 


To determine whether a point P(x, y, z) (such as a control point) is 
within the boundary of the plane formed by the points 1-2-3-4, the two 
vectors and $2 formed by points 1-P and 1-2, respectively, can be 
used. When projected onto the X— Y plane, these two vectors will be 
denoted by and Q^- It follows that 


Qi X Q’ = 


1 

X - X, 


x^ - x^ 


j 

y - 


Y2 - y 


k [(x - Xj^) (y 2 - y^^) - (X 2 - x^) (y - y^)j 


(B.8) 


It follows that if x > 0, the point P is on the flap-side of the line 
connecting points 1 and 2. In case is a zero vector (i.e., points 1 and 
2 coincide), comparison of y-coordinates of points 1 and P will be sufficient 
to determine whether P is inboard or outboard of point 1. Similarly, the 
two vectors Pj^ and P^ formed by points 4-P and 4-3, respectively, and projected 
onto the X-Y plane can be used to determine the location of P relative to 
the boundary of the line 4-3; 
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j 

y - 74 
73 - 74 

= S^[(x - X^)(y 3 - 7^) - (X 3 - x^)(y - y^)] (B.9) 

If X < 0, the point P is again on the flap-side of the line connecting 
points 4 and 3. Again, if points 4 and 3 coincide, the y-coordinates of 
point P and 4 may be compared to determine the location of P relative point 4. 

B.2 Conical Camber 

Leading-edge plane flaps, described in Section B.l, are simple devices 
frequently used to simulate aerodynamically the effect of more complicated 
camber shapes, such as conical camber. However, many test models have been 
incorporated with leading-edge conical camber of different sizes. To illus- 
trate the derivation of analytical expressions describing the geometry of 
such conical cambers, the one used in reference 7 on a delta wing will be 
used as an example. 

As shown in figure B.2, Q^Q 2 is s circular arc with circle center at 
O'. Since the point is on the circular arc, it follows that (see fig. 

B.2(a)) 

(0.195a)^ + (B.IO) 

0.102a + z| = R (B.ll) 

Substituting equation (B.ll) into equation (B.IO) results in 

(0.195a)^ + (R-0.102a)^ = R^ 
which can be solved for R; 

R = 0.237397 a (B.12) 

Based on the planform coordinate system as shown in figure B.2(b), the 
center of circle O' is at (x, y, -R) , where 



k 

0 

0 
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y = 0.231303 (1.6+x) 


(B.13) 


R = 0.237397 a = 0.237397 x (1.6 +x) 

= 0.06821207 (1.6+x) (B.14) 

Hence, the equation of circular arc based on the planform coordinate system 
becomes 

(y - y)^ + (z + R)^ = R^ , y>^y , z < 0 (B.15) 

3z 

To find replace z in equation (B.15) by z^ and differentiate the 

equation with respect to x. It is obtained that 

2(y - y) (- g) + 2(i^ + R) (-^ + g, . 2 ,|S 

but 

^ = 0.231303 
3x 


= 0.06821207 
3x 


Hence , 

3 z 


3x z + R 
c 


[ 0.231303 (y - y) - 0.068212 (z^)], y ^ y (B.16) 


where z is given by 


z^ = -R + J~R^ - (y - y)^ , y y 


(B.17) 


3z 

To find — — , equation (B.15) is differentiated with respect to y: 
3y 


3z 
c 

3y 


y - y 

Z + R 
c 


y >■ y 


(B.18) 
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APPENDIX C 


A Numerical Method for Calculating 


3z 
c 

3y 


When a general camber surface is described by a finite number of ordinates 

3z 

c 

at some spanwise stations, a numerical method is needed to calculate -r — and 

— oX 

9z 


at control points and integration (pressure) points. In the present 
3z 


3y 

u ^ 

program, -T — is calculated through cubic spline interpolation. In the 
ox 


following, a trigonometric interpolation formula will be used to calculate 
3z 

The method used is similar to that used in determining the distribution 
3y 

of streamwise vortex density (ref. 12). 


It is well-known that if at a given point, slopes in two directions are 
known, then the slope in other directions can be determined in terms of them. 
Referring to figure C.l, assume that the slopes in the X- and directions 
are known. To calculate the slope in the y-direction, the following 


coordinate transformation formulas are needed. 
X* = X cos A - y sin A 
y'=xsinA+ycosA 


If z^ (x,y) is the dimensional 


camber function, then 


9z 

c 

9z 

_ c 

3x' 

3z 

+ ^ 

3y' 

9y 

9x^ 

3y 

3y' 

3y 

9z 

c 

9z 

c 

3x' 

3z 

, 9. 


9x 

9x* 

9x 

3y' 

3x 


-sin 


3z 

A T + cos A 

3x 



3z 


3z 


c c 

cos A -r — r + sin A -r — r 
3x 3y 


(C.l) 

(C.2) 


(C.3) 


(C.4) 


Substituting — obtained from equation (C.4) into equation (C.3) results in: 


3z 3z , 3z 

c , c , 1 c 

-r = -tan A 1 r “r — r 

3y 3x cos A 3y 


(C.5) 


Equation (C.5) represents the basic formula used in the present method. To 
3z 

determine -r — r*z -ordinates in the y '-direction are used. Although cubic 


3y 


3z 


spline interpolation can be used to find "g^rr was found to be more 


accurate to use the trigonometric interpolation formula (ref. 12): 
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(G.6) 


3z dz „ 

c _ c I 

8yl d9 . b' sin 8. 
'^11 1 


dz 


de . 

X 


M-1, 

z (-1) 
k=l k 


, , . sin 6 . , cos 0 . 

k+x X , 1 X 

cos 0, - cos 0. 2 c. sin 0. 

k XXX 


. (-1) sxn 0. 

,1 X 

2c 1 - cos 0 . 


“ 2 ^ 


Cm 


(- 1 ) 


M+i 


sin 0 . 
3^ 

1 + cos 0 . 


where : 

y' = ^ (1 - cos 9) 


0 . 
X 


ir 

M 


1 ,.. 


,M-1. 


kti 

M 


(C.7) 


(C.8) 

(C.9) 


where the prime on the summation sign implies that the term with k = i is to 


be omitted. In equation (C.7), z and z are the camber ordinates at 
inboard and outboard endpoints, respectively, and M is the total number of 
points used in the interpolation. 

To indicate the accuracy of the interpolation method, an arrow wing 
with the following camber surface is chosen (See figure C.2): 

(x, y) = (0.14 - 0.22y)J^ + (-0.1835 + 0.165y)J^ + 0.0785x 


+ 0.015 (y - 1.0) 
y = y/b/2, x 


(x - (y))/c(y), z^ = z^/c(y) 


It follows that 
3z 3z 


3y 


C 

3y 3y 


(C.IO) 

(C.ll) 


Iz 

3y 


[c(y) + z^ 


3c 


(C.12) 
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For the configuration shown in figure C-2 


X = (X - y tan A^)/c(y) = 


where ; 


c(y) = (1 - y) 


3c(y) 

3y 


-c = -1.95 
r 


From equation (C.IO), 

9 z _ — „ 

= -0.22 + 0.165X + 0.015 + [(0.42 - 0.66y)x 


+ (-0.367 + 0.33y)x + 0.0785] 

Using equation (C.13), the derivative can be obtained: 

dy 


3x _ X - 3 
3y ~ c^ (1 - y)2 


Substituting equations (C.14)-(C.17) into (C.12) results in 

3z ? r _ 

— - = — 1.95 (1 - v) — - - 1.9“= 

3y b L ^ 3y c 


(y)] 


(C.13) 


(C.14) 

(C.15) 


(C.16) 


(C.17) 


(C.18) 


with 


3z 3 9 

-=^ = -0.22x + 0.165x + 0.015 + -s 

3y 1.95 (1 - y) 


(x - 3 ) (0.42 - 0.66y) - 2 


(x - 3) (-0.367 + 0.33y) - x - 3 

+ 1.95 (1 - y)2 ^ 1.95 (1 - y)'2 


(C.19) 

3z 


Equations (C.IO), (C.18) and (C.19) describe the exact distribution of 

dy 

for the assumed arrow wing geometry. 

Numerical results at three spanwise stations are compared with exact 
values in table C.l. The results along the leading edge are compared in 
table C.2. From both tables, it can be seen that the present numerical 
method is quite accurate. 
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3z 

Q 

Table C .1 Comparison of Calculated - 5 — with Exact Results 

dy 

at Three Spanwise Stations 


y 

X 

X 

3zc/ 3y 
Calculated 

Exact 

.172570 

.588089 

.043620 

-.135850 

-.135850 


.783837 

. 164940 

-.054444 

-.054436 


1.092941 

.356414 

.023043 

.023043 


1.448850 

.577098 

.036728 

.367278 


1.781067 

.782997 

-.021320 

-.021320 


2.023794 

.933434 

-.105749 

-.105749 


2.128961 

.998613 

-.153111 

-.153111 

.571157 

1.749949 

.043620 

-.171606 

-.171606 


1.851402 

.164940 

-.116446 

-.116445 


2.011605 

.356414 

-.042382 

-.042382 


2.196066 

.577098 

.017433 

.017433 


2.368248 

.782997 

.042307 

.042307 


2.494050 

.933434 

.037956 

.037955 


2.548556 

.998613 

.029470 

.029470 

.979746 

2.940962 

.043620 

-.208260 

-.208256 


2.945754 

.164940 

-.180003 

-.179998 


2.953320 

.356414 

-.109448 

-.109446 


2.962032 

.577098 

-.002346 

-.002347 


2.970163 

.782997 

.107531 

.107529 


2.976105 

.933434 

.185266 

.185260 


2.978679 

.998613 

.216633 

.216628 
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with Exact Results 


Table c . 2 


Comparison of Calculated 


Sz 
c 

9y 


Along the Leading Edge x=0 


y 

y 

X 

3z / ^ 
c 

Calculated 

Exact 

.022115 

.020254 

.060761 

-.163188 

-.163188 

.086668 

.079373 

.238120 

-.166355 

-.166355 

.188429 

.172570 

.517709 

-.171348 

-.171349 

.319154 

.292292 

.876877 

-.177763 

-.177764 

.468253 

.428843 

1.286528 

-.185079 

-.185080 

.623647 

.571157 

1.713472 

-.192703 

-.192706 

.772746 

. 707708 

2.123122 

-.200019 

-.200024 

.903471 

.827430 

2.482291 

-.206433 

-.206444 

1.005253 

.920627 

2.761880 

-.211427 

-.211452 

1.069785 

.979746 

2.939239 

-.214594 

-.214698 
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Figure 1 Definition of Variables Defining a Nonplanar Wing Surface. 
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Leading; Edge 
IV 


Figure 6 Measured Variation of Circumferential Velocity 

in Four Quadrants (ref. 10) at Angle of Attack of 
14.9 deg. for A=1.0 Delta Wing, 






Figure 7 Geometry and Flow Field for Defining Vortex 
Center Location , 



Figure 8 Velocity Components along the Camber Surface 


with a 


0 . 












95.964 

(37.781) 



739 . 54 ^ 

110.754- 

(43.604) 


Figure 12 Geometry for a Supersonic Cruise 
Dimensions in cm. (in.) 


FLAT DIHEDRAL 
MODEL MODEL 






ration Defined in Reference 15 









Figure 15 Effect of Vortex Lift Orientation on Predicted Aerodynamic 
Characteristics of A=1.147 Delta Wing with Conical Camber 
at M=0. 




Figure 16 Vortex Locations on A=1.07 Delta Wing 
= 12.29 deg. 


at a 



L 






18 , ( 



Figure 19 Geometry of a 74°-Delta Wing with Leading 
Edge Vortex Flap Defined in Reference 8, 















O Data (ref.23) 


Present Theory with Correction, Vortex Flow 

Present Theory without Correction, Vortex Flow 

Present Theory, Attached Flow 

Original Suction Analogy (ref. 22) 




Figure 26 Aerodynamic Characteristics of a Plane A=1.333 
Delta Wing at Ii=1.8. 
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Figure 27 Angle of Attack Correction for Supersonic Vortex 
Separated Flow at a=10 deg. 



Figure 28 Angle of Attack Correction for Supersonic Vortex 
Separated Flow at Different a and M . 
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